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ABSTRACT 

We study the X-ray emission from five symbiotic stars observed with Suzaku. These objects 
were selected for deeper observations with Suzaku after their first detection with ROSAT and 
Swift. We found that the X-ray spectra can be adequately fit with absorbed optically thin thermal 
plasma models, either single or multi-temperature. Such a model is compatible with the X-ray 
emission being originated in the innermost region of the accretion disk, i.e. a boundary layer. 
Based on the large flickering amplitude (only detected in 4 Dra), the high plasma temperature 
and previous measurements of UV variability and luminosity, we conclude that all five sources 
are accretion-powered through predominantly opticall thick boundary layer. 

Given the time lapse between previous and these observations, we were able to study the 
long term variability of their X-ray emission and found that the intrinsic X-ray flux and the 
intervening absorption column can vary by factors of three or more. However, it is still elusive 
the location of the absorber and how the changes in the accretion rate and absorption are related. 

Subject headings: binaries: symbiotic. X-rays, individuals: CD -28 3719, EG And, Hen 3-461, Hen 3-1591, 
4 Dra 


1. Introduction 

When observed at optical wavelengths, symbi¬ 
otic stars (SS) show a composite spectrum which 
suggests that they are binary systems. A hot, 
compact component (usually a white dwarf, WD) 
contributes to the blue-UV region of the spec¬ 
trum, while a cool red giant dominates the spec¬ 
trum at longer wavelengths. Observations at other 
wavelengths reveal a very complex and rich sce- 
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nario for these systems. Optical, infrared and UV 
spectral regions are rich in emission lines from 
forbidden and permitted transitions, which arise 
mainly from photoionization and recombination 
of the n ebular plasma heated by the the hot com- 
Done nt dKenvon 2009) Radio observations reveal 
jets ( BrocksonneLaL- 2004ICrocker et iiil~ll2001 


Kellogg et al. 2007t) and the rmal emission from 


the i onized red-giant win d (|Seaauist fc Tavlor 


1990l : Seaauist et al.l[l9^ . Symbiotic stars can 
even produce y-rays during n ova-type outbursts 
(e.g. V407 Cyg; V745 Sco lAckermann et al 

20141) . 


Symbiotics are now recognized as a popula¬ 
tion of X-ray sources. From the ~ 220 systems 
known, 45 have been detected at X-ray wave¬ 
lengths, most of them with emission in the 0.3 
to 10 keV range. A few, h owever, were detecte d 


at energies up to 100 keV ( Kennea et al. 20091) . 


The X-ray spectra from sy mbiotic stars c an be 
divided into five categories ( Luna et al. ( 2013h 


hereafter Paper 1), according to the nature of 
their emission (thermal or non-ther malj and the 
maximnm energy of their spectra. iLuna et al. 


1 












































( 2013 ) and Nunez at al. ( 20141) studied and clas¬ 
sified the first X-ray detections of many symbiotics 
with Swift, XWIM-Newton and/or Chandra. They 
found that thermal X-ray emission (optically thin 
or blackbody-type) arises from the accretion-disk 
boundary layer, a colliding-wind region or quasi¬ 
stable nuclear burning on the WD surface. 

Sensitive, broadband X-ray satellites such as 
Suzaku have also played a significant role in 
observing symbiotics, especially those with X- 
ray emissi on above 10 keV (T CrB, CH Cyg . 


ray emissi on above lU keV (1 (Jrtj, CM Uyg 
V648 Car; Luna et al. Il2005 iMiikai et al. 11200^ 
Kennea et al. Il2009[) . Here we present the analysis 


of Suzaku/NlS observations of a set of five symbi¬ 
otics previously known to be X-ray sources. This 
sample consist of CD-28 3719, EG And, Hen 3- 
461, Hen 3-1591 and 4 Dra. In Section[2]we detail 
the data reduction and analysis, while in Section 
[3] we present the results. Concluding remarks are 
presented in Section H) 


2. Observations and data reduction 


The X-ray Imaging Spectrometer (XIS) on 
board Suzaku observed five symbiotic stars. De¬ 
tails of each observation are presented in Table [I] 
All sources were observed with the XISO, XIS I, 
and XIS3 detectors ( the XIS2 is no t operative 
since November 2006 (jKovama 1120071 1. which are 
sensitive in the 0.2-12.0 keV range. All sources 
were too faint to be detected with the Hard X-ray 
detector (HXD). We reprocessed all data using 
the aepipeline script and obtained event files 
with the processing version 2.5.16.29 (2014-07-01) 
applied. 


Source spectra and light curves were extracted 
from circular regions centered on the source SIM- 
baeQ coordinates. The recommended radius of 
the extraction regiorU is 260" (this size encircles 
99% of the point source flux), however given the 
source and background brightness, we were able 
to use this size only in the case of 4 Dra. Compar¬ 
ing the source spectrum with that from the back¬ 
ground, we found that the optimal radius for the 
source region, which maximizes the signal-to-noise 
ratio, was of 120" in the case of Hen 3-1591, CD - 
28 3719 and EG And and 60" for Hen 3-461. Back- 


^ http: / / simbad.u-strasbg.fr / simbad/sim-fid 
^http://heasarc.gsfc.nasa.gov/docs/suzaku/analysls/abc/ 


ground spectra and light curves were extracted 
from annulus regions centered on the source (with 
inner and outer radius of 200" and 400", respec¬ 
tively) in the case of Hen 3-1591, CD -28 3719, 
Hen 3-461 and 4 Dra while a circular region with 
160" radius was used for EG And because the lo¬ 
cation of the source on the chip did not allowed to 
select an annulus region for the background. The 
response matrices were created using the rmfgen 
and arfgen scripts. We then fit the binned spec¬ 
tra (grouping by a minimum of 16 to 20 counts 
per bin) using XSPECH and as statistic to se¬ 
lect the best fit model. 


Hen 3-1591 was observed serendipitously with 
ASCA ( Tanaka et al. 1994[) on 1999 September 
22 during an observation of the supernova remnant 
(SNR), G5.2-2.6. The Solid-state Imaging Spec¬ 
trometer (SIS) was operated in I-CCD mode for 
this observation (of the 4 chips available), which 
put Hen 3-I59I outside the operational SIS field of 
view. On the other hand. Hen 3-I59I was securely 
in the field of view of the Gas Imaging Spectrom¬ 
eter (GIS) instrument, which has two units (GIS2 
and GIS3). We have selected intervals when the 
satellite was outside the South Atlantic Anomaly 
(SAA), the attitude control was stable and within 
0.02 degrees of the target, and the line of sight was 
greater than 5 degrees above the Earth limb. We 
have also applied a standard selection expression 
combining monitor count rates and geomagnetic 
cutoff rigidity (COR) to exclude time intervals of 
high-particle background, obtaining 15 ks of good 
on-source data. 


Eor the ASCA observation on Hen 3-1591 we 
have extracted the source spectrum from a 6' ra¬ 
dius circular extraction region centered on the 
source, and background spectrum from a 6' radius 
region centered at (271.6619, -25.8192) away from 
the SNR and other obvious sources. Hen 3-1591 is 
detected at a net count rate of 0.016 c/s/GIS. We 
use the v4.0 rmf downloaded from the CALDB and 
have generated an arf file for each unit of the GIS. 
We have then combined the GIS2 and GIS3 spec¬ 
tra and responses, ignored the data outside the 
well-calibrated range of 0.7-10 keV, and binned 
the data by a factor of 32, leaving 26 channels. 

We evaluate whether any of the three spectral 
models, i.e. an absorbed single-temperature op- 


® http: //heasarc.gsfc.nasa.gov/docs / xanadu/xspec / 


2 



























Table 1 

Observations Log 


Source 

Date 

Obsid 

Exp. Time [ks] 

4 Dra. 

2010-04-18 

405035010 

42 


2011-11-09 

406041010 

42 

CD -28 3719 . . 

2013-10-12 

408032010 

14 

EG And .... 

2011-02-05 

405034010 

100 

Hen 3-461 . . . 

2012-12-17 

407007010 

46 

Hen 3-1591 . . 

2012-10-03 

407042010 

51 

- ASCA .... 

1999-09-22 

57055000 

15 


tically thin thermal plasma, an absorbed multi¬ 
temperature cooling flow or an absorbed non- 
thermal plasma properly fits the data by their 
xl- Those best-fit models are described in the fol¬ 
lowing sections while the resulting parameters are 
listed in Tabled With the selected model, we use 
the unbinned data and C-statistic ICash (197^ 
to calculate the uncertainties in the parameters of 
the model and the flux. All errors in the fit pa¬ 
rameters were estimated at their 90% confidence 
(see Tabled]). 


3. Results 

3.1. Spectral analysis 

The five sources observed with Suzaku were 
modeled with absorbed optically thin thermal 
emission, which in the case of Hen 3-1591, sug¬ 
gest that the compact object is a white dwarf 
instead of a neutron star as previously thought. 
All sources show strong changes in their X-ray 
flux and the hardness in long time scales, from a 
few months to years, which could be related with 
changes in the accretion rate and/or absorption. 


3.1.1. 4 Dra 


This source was classified as a triple sys- 
tem cons i sting of 4 Dra(A) -|- CQ Dra(Bab) 
( Reimers li^Sll . However, X-ray data o btained 
with ROSAlQlead Wheatley et al. I ( 2003ll to sug¬ 
gest that the X-ray emission is consistent with 
the presence of a WD accreting from the wind 
of the r ed gia nt and the UV data from lUE led 
Skopali (2005) to the same conclusion after fit¬ 
ting the SED. More evidence against the triple 


system nature of 4 Dra comes from the analy¬ 
sis of the broad wings superimpo sed to narrow 
emiss ion lines in FUSE spectra by ( Froning et al.l 
I 2 OI 2 II . similar to other FUSE spectra of confirmed 
symbiotic stars. This source is one of the two 
system with H ipparcos distance in our sample, 
fi=190±17 pc ( van Leuween I [2007 1. Its optical 
brightness presents irre gular variations about 0.1 
mag in V ( Eggen I[l967ll and have radio flux varia¬ 
tions on time scales of weeks to months an d maybe 
shorter than hours to days (iBrown 111987^ . Based 
on these previous studies, in order to model the 
Suzaku spectra, we consider this source to be a 
symbiotic system with an accreting white dwarf. 


We analyzed the two X-ray observations, Ob- 
sID 405035010 and 406041010 (hereafter 4050 and 
4060, respectively. See Table [T]). We obtained ac¬ 
ceptable fits with two spectral models in two dif¬ 
ferent scenarios. In the first scenario, we forced the 
temperature of the optically thin thermal emission 
to have a unique value for both observations. This 
is a valid assumption if the X-ray emitting plasma 
arise in the post-shock region of the accretion disk 
boundary layer and its temperature is set by the 
WD mass, which should not change between ob¬ 
servations. In the second scenario, we modeled 
both observations independently. 


The first spectral model consists on a sin¬ 
gle temperature APECi plasma model with a 
reduced metalicity, observed through a simple 
and partial covering absorber. For the first sce¬ 
nario the plasma temperature and abundances 
are A:T=6.0±0.3 keV and abundance=0.18±0.02; 
the absorption column of the simple absorber is 
^ff °22 =0-43±0.03 and Af^°f^=0.94±0.06; while for 


“^http: //www.xray.mpe.mpg.de/cgi-bin/rosat/rosat-survey 


® https: //heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XSmodelApec.html 
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the partial covering absorber we have A^^°|2=l-26±0.5, 3.1.2. CD -28 3719 


Ff5O=i30±2x 10-13 ergs cm-^ 
ergs cm”^ s 


H,22- 

-^_h° 22=2.61±0.22; covering fractions of CF(4050)=0.43±0.03 
and CF(4060)=0.64±0.03; unabsorbed fluxes 

s-i,Ff30=697±2x 10-13 

1 and Xy=l-06. 

If we use the model above to fit both observa¬ 
tions independently, we obtain A:T(4050)=4.5±0.1 
keV and fcr(4060)=6.4±0.1 keV; abundances of 
0.17 (4050) and 0.19 (4060); iV^°^o=0.46±0.03 and 
while for the partial-covering 


lV^o^§=0.88±0.07 
absorber we have A'^^3|2=2.5±0.5, Nfj 22 =‘^-^^ 0 . 2 ; 
covering fractions of CF(4050)=0.38±0.04 and 


In X-rays, CD -28 3719 was detected for the first 
time in a short pointing observation with Swift and 
its spectrum classified as 5-type, i.e. optically thin 
thermal, hard X-ray emission observed through 
strong absorption (Paper I). Our Suzaku obser¬ 
vation aimed to obtain a spectrum with a higher 
signal-to-noise ratio in order to improve the basic 
spectral modeling. During the observation there 
were some problems in the acquisition of the XISO 
chip data which did not go back to the 5x5 editing 
mode after dark frame dump during a South At- 


CF(4060)=0.66±0.04; unabsorbed fluxes F^30_236±2xl(laMic Anomaly passage and the data during these 

segments were corrupted. For this reason, we only 
analyzed the data of the XISl and XIS3 chips. 


ergs cm 

,-i 


-2 „-l 


F^030=564±2x10-i3 


ergs cm 


-2 


s-i and x2(4050)=1.02, (4060)=1.07. 


As an alternative, the second model consists 
on a multi-temperature cooling flow model ob¬ 
served through a simple and partial covering ab¬ 
sorber. In the case that the X-ray emission from 
both observations originates from a plasma with 
the same temperature and abundances, we ob¬ 
tain fcT=I0.8tQ5 abundance=0.2I±0.0I; 

the absorption column of the simple absorber is 
A4050=o.54±0.07and A^^°3°=1.21±0.06; while for 
the partial-covering absorber we have 
-^ff°22=3.5±0.3; covering fractions of CF(4050)=0.391 q Jg 
and CF(4060)=0.60±0.03; unabsorbed fluxes 
Ff30=145±2x 10-33 ergs cm-^ s-^, Ff30=723±2x 10-33 
ergs cm-^ s-3 and x^=1.05. 

If we model both observations independently 
with an absorbed cooling flow model, we obtain 
fcr(4050)=7.3tj;^ keV and fcT(4060)=12.9;)!j 
keV; abundances 0.18±0.03 (4050) and 0.22±0.02 
(4060); iV^030^o.64±0.05 and 7V^030=1.13±0.07; 
while for the partial-covering absorber we have 
22=3.2±0.9, iV^3f2=2-8±0.3; covering frac¬ 
tions of CF(4050)=0.42±0.04 and CF(4060)=0.61±0.04; 
unabsorbed fluxes F^33=175±2x 10-33 ergs cm-^ 
s-3, F^33=673±2xl0-33 ergs cm-^ s-3 and 
x2(4050)=1.01 x^(4060)=1.07. 

Although in terms of xt s-n absorbed non- 
thermal plasma plus a gaussian emission line 
(wabsx (power+gauss)) fit the observed spec¬ 
trum, the line centroid is at ~6.67 keV, consistent 
with Fe XXV transitions from a thermal plasma. 

We thus opt for a thermal origin for the observed 
X-ray emission. 


The best fit model for this source consist on an 
absorbed cooling flow model with variable abun¬ 
dance wabsxmkcflow. The absorption column 
Nh ,22 is high, 12 . 9^2 g! maximum temperature 
kTjnax is 8 I 3 keV and the abundance is 0 . 44 to;i 7 . 
The unabsorbed flux (in the 0.3-10 keV range) Fx 
is 23 I 2 XIO -33 erg cm-^ s-3. The luminosity Lx, 
at a distance of 1 kpc, is 28^2 x 10^^ ergs s-3 (d /1 
kpc)^. The results from this fit are commensu¬ 
rate with the results obtained from the Swift data 
analyzed in Paper I. 


3.1.3. EG And 


In the V band, EG And is one of the bright- 
est symbiotic system s. Its distance is 5I2±168 pc 


(Ivan Leuween 1200711 . The periodic photometric 


modulation indicates an orbital perio d of 470 days 


with an inclination of 82/^® g degrees ( Kolb et al. 


20041 1. making EG And an eclipsing symbiotic bi¬ 
nary. Our Suzaku observation took place at the 
orbital phase (5= 0.93 (using the ephemeris from 


Kolb et al. 1120041) , with the WD going behind the 


red giant wind or (p=0.17 is we use the orbital pe¬ 
riod of 481 days and the ephemeris from Vogel et 
al. (1991), i.e. the WD coming out of the partial 
eclipse. 

After visual inspection, we verified that no 
obvious features were present in the spectrum 
and group the channels to have a minimum 
of 20 counts per bin, which allow us to use 
the x^ statistics during the fit. The best fit 
model consists of an absorbed single plasma tem¬ 
perature wabsx ape c, with absorption column 
Nh.22=0-0^-o'^j and temperature kT= keV. 
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The unabsorbed flux is Fx=l-3±0.2xl0 ergs 
s“^ cm“^ and the luminosity Ljf=0.4x10^^ ergs 


3 . 1 . 4 . Hen 3-461 

Hen 3-461 was discovered in X-rays with Swift 
during a short pointing observation (Paper I). The 
high temperature and absorption obtained from 
modeling those data, the hardness ratio (defined 
as the ratio of count rates in the 2.4-10/0.3-10 
keV ranges) and the presence of significant flick¬ 
ering in the UV, suggested that the X-ray emis¬ 
sion originates in the boundary layer of the ac¬ 
cretion disk and lead us the authors classify it 
as a 5-type source. Almost two and half years 
later, Suzaku observed Hen 3-461 and only de¬ 
tected 285 source photons, i.e. 8 cr detection af¬ 
ter background substruction. Given the low num¬ 
ber of photons detected, we performed a crude 
spectral modeling (x^=1.45). We noticed that a 
weak soft component is detected at a low signif¬ 
icance level ( 2 cr) and thus we decided to include 
a second component in our spectral model, being 
aware of the similarity with the well-known, two- 
spectral components /3/5-type emission observed 
in few white dwarf symbiotics (e.g. NQ Gem, 
V347 Nor, see Paper 1). Our best spectral model 
consist of two optically thin thermal components 
(fcri=0.07±-k0.02 keV and kT 2 =^tl keV) and a 
gaussian emission line at 6.4 keV. The hardest 
component is modified by a full (iVij^ 22 = 3 i 2 ) 
a partial covering absorber ‘hov¬ 

ering fraction=0.97±0.02). The unabsorbed X- 
ray flux is Fx=26±5xl0~^^ ergs s“^ cm“^ and 
assuming a distance of 1 kpc, the luminosity is 
Lx=31±lxl0^^ ergs s“^. 


3.1.5. Hen 3-1591 


On ly two objects were classified bv lMiierset et al 
( I 997 II as 7 -type symbiotics based on their emis¬ 
sion detected with ROSAT; GX 1+4 and Hen 3- 
1591. These o bservations had limite d signal-to- 
noise ratio but iMiierset et al. (jlQQTI l speculated 
that these sources harbored neutron stars as accre- 
tors. Lately, these type of symbiot ics were named 
as ’’sy mbiotic X-ray binaries” by iMasetti et al. 
(I2nn7ll . 


We have found that Hen 3-1591 was serendipi- 
tously observed with ASGA/GIS (but outside the 


SIS field of view) in 1999, with an exposure time 
of 16 ksec. The spectrum is not of high enough 
statistical quality (see Fig. [2|) to draw definitive 
conclusions (only ^200 photons were detected), 
but some inferences can be drawn. A power-law 
fit prefers a photon index near 1.9, with excess 
counts around 6.5 keV. Adding a Gaussian, the 
line centroid was found to be near 6.6 keV with 
an equivalent width well in excess of 1 keV. We 
interpret this as due to a combination of Fe H-like, 
He-like, and fluorescent lines, i.e., the X-ray emis¬ 
sion should be modeled as optically thin thermal 
emission with reflection adding the 6.4 keV line. In 
fact, a cooling flow plus Gaussian model fit gives 
a good description of the observed spectrum. The 
maximum temperature is around 14 keV, iron is 
strongly overabundant (~ twice the solar value), 
and the equivalent width of the fluorescent line 
600 eV) also requires an overabundance of iron. 
Taken at face value, this suggests that Hen 3- 
1591 hosts a white dwarf accreting matter over¬ 
abundant in Fe. Our Suzaku spectrum aimed to 
identify, in a higher signal-to-noise spectrum, the 
best spectral model. Given that ROSAT data sug¬ 
gested a non-thermal origin for the X-ray emission 
while ASCA data suggested that the X-ray emis¬ 
sion is due to a multi-temperature optically thin 
thermal plasma, we try these two models in the 
Suzaku spectrum. First we try two thermal mod¬ 
els, which show the same statistic (xJ=l-07) and 
similar model parameters. We prefer the multi¬ 
temperature model, wabsxmkcf low, which yielded 
an absorption column of A^//, 22 = 0 . 1 l//o;o 6 and a 
maximum temperature kTmax=7-2 keV. The un¬ 
absorbed flux is Fx=7.4±0.3x 10“^^ ergs cm“^ 
s“^, and Lx at a distance of 1 kpc is 8.9+0.4x10^^ 
ergs s“^. 

On the other hand, an absorbed non-thermal 
model, wabsxpowerlaw yields Afi/_ 22 = 0 . 21 ± 0.06 
and the photon index is 2.40+0.19. The ab¬ 
sorbed flux Fx=10.7+0.5x10“^^ ergs cm“^ s“^, 
and the luminosity Lx=12.9+0.6xl0^^ ergs s“^ 
(d/1 kpc)^ x^=l.ll. 

3.2. Timing analysis 

Significant stocasthic variability of the the X- 
ray flux in short (minutes to hours) time scales is 
a hallmark of the presence of accretion disk in bi¬ 
nary systems. Periodic modulation, on the other 
hand, would indicate that the accretion is chan- 


5 

















10-3 


10 ^ 





'tfjp'yp 


pf 


Energy (keV) 




0.5 1 2 5 10 

Energy (keV) 




Energy (keV) 


Fig. 1.— Suzaku/XlS) spectra of (a) CD -28 3719, (6) Hen 3-1591, (c) Hen 3-461, {d) EG And and (e) 4 Dra 
ObsID 4050 and (/) 4 Dra ObsID 4060. The full line shows the best-fit model described in Section [3l Red 
(XISl), black (XIS0-h3) 
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Table 2 

Spectral Models 


Object 


Model 


Count rate Nh ,22 

[10-2 counts s“^] 


kpa p^b 

[keV] 


1st scenario 


(see Sect. 13.1.Ill 


4 Dra(4050) ... 

wabs X pcfabs xvapec 

38±1 

0.43±0.03 

6.0±0.3 

136±2 

6±1 

1.06 




1.26±0.50, CF=0.43±0.03 





4 Dra(4060) ... 

wabs X pcf abs X vapec 

110±1 

0.94±0.06 

6.0±0.3 

697±2 

30±1 

1.06 




2.61±0.22, CF=0.64±0.03 





4 Dra(4050) ... 

wabs X pcfabs X mkcflow 


0.54±0.07 

I0.8l“;« 

145±2 

6±1 

1.05 




1.5±0.5. CF=0.39l“;“ 





4 Dra(4060) ... 

wabs Xpcf abs xmkcf low 


1.21±0.06 

io.8l«;« 

723±2 

31±1 

1.05 




3.5±0.3, CF=0.60±0.03 





2nd scenario 








(see Sect. |3.1.ip 








4 Dra(4050) 

wabs X pcf abs X vapec 


0.46±0.03 

4.5±0.1 

136±2 

6±1 

1.02 




2.5±0.5, CF=0.38±0.04 





4 Dra(4060) 

wabs xpcf abs xvapec 


0.88±0.07 

6.4±0.1 

564±2 

24±1 

1.07 




2.3±0.2, CF=0.66±0.04 





4 Dra(4050) ... 

wabs xpcf abs xmkcf low 


0.64±0.05 

7±1 

175±2 

8±1 

1.01 




3.2±0.9, CF=0.42±0.04 





4 Dra(4060) ... 

wabs xpcf abs xmkcf low 


1.13±0.07 

12±1 

673±2 

29±1 

1.07 




2.8±0.3, CF=0.61±0.04 





CD -28 3719 ... 

wabs X apec 

0.9±0.1 


5±2 

17±1 

20±1 

1.17 


wabs xmkcf low 


13±3 

o -|-6 

®-3 

23±1 

28±2 

1.21 

EG And . 

wabs X apec 

0.4±0.1 

u . uy_Q 07 

7+^ 

‘ -2 

1.3±0.2 

0.4±0.1 

0.9 

Hen 3-461 . 

apec+wabs xpcf abs X (apec+gauss) 

0.2±0.1 

q -|-5 

'^-2 

q -|-5 

°-3 

26±5 

31±6 

1.45 




75±25; CF=0.97±0.02 





Hen 3-1591 .... 

wabs X apec 

3.7±0.1 

0.05±0.04 

3 . 06 ; s :48 

a ^-|-0 . 3 

6 - 6 _ o .3 

7.9±0.4 

1.07 


wabs xmkcf low 


0 lX + 0.05 

7 - 1-3 
‘ -2 

7.4±0.3 

8.9±0.4 

1.07 


®^Indicates the value of the maximum temperature in the case of the cooling flow model mkcflow, kTm, 
^Unabsorbed X-ray flux, in units of 10”^^ erg s~^ cm~^ in the 0.3-10.0 keV energy range. 
‘^Unabsorbed X-ray luminosity, in units of 10^^ erg s“^ in the 0.3-10.0 keV energy range. 

‘^Absorption column and covering fraction of the partial absorber model pcfabs. 


7 
















neled by an strong m agnetic field (e.g. Z And, 


Sokoloski et al. ( 2006l l'l. 


We searched for periodicity using the Lomb- 
Scargle algorithm on the light curves binned at 
16 s. No periods were found for the sources in our 
sample, with a false alarm probability >0.3. We 
used the ratio between the observed (s) and ex¬ 
pected (Sexp) variances to account for stocasthic 
variations in the light curves. We binned the 
lightcurves at integer numbers of the XIS readout 
time: 16, 160, 1600, and 3680 seconds and calcu¬ 
lated the ratio sjsexp (see Table E]). Only 4 Dra 
shows significant variability, with amplitudes as 
high as 40%, the strongest flickering is present in 
hours scales but there is also intense flickering in 


minutes time scale. 


terminations (jWheatlev et al 


consistent with previous de- 

200il . 


4. Discussion and Conclusions 

We analyzed the deep, broad-band Suzaku 
observations of five symbiotic stars, previously 
known to be X-ray sources. The spectra obtained 
with Suzaku allow to particularly unveil the origin 
of the X-ray emission in symbiotics when com¬ 
pared with data obtained with ROSAT, limited 
by energy coverage, or Swift, limited by contin¬ 
uous observing time. The high temperatures of 
the plasma (see Tabled strongly suggest that the 
emission originates in an accretion disk bound¬ 
ary layer (5-type) instead of a colliding wind re¬ 
gion (/3-type). Of the five sources observed with 
Suzaku, only CD -28 3719 retains its previous 
classification as 5-type source (see Paper I) de¬ 
rived from Swift observations. The proposed WD 
nature of the compact object in Hen 3-1591 and 
the temperature of the X-ray emitting plasma 
suggest that it should now be considered a 5- 
type source instead of y-type source (those sym- 
bio tics with neutron st ars as compact objects; 
see iMiierset et al. Il997r i . The likely presence of 


a new soft component in the X-ray spectrum of 
Hen 3-461 encourage to propose a /3/5-type clas¬ 
sification. No X-ray spectral type was proposed 
for 4 Dra and given the results obtained from the 
spectral fit, we advocate for a 5-type categoriza¬ 
tion. Finally, the broadband energy coverage now 
shows that EG And should be considered a 5-type 
instead of /3-tvpe sour ce as originally proposed by 
Miierset et al. I (1997]). While the 5 components 


in the iKennea et al. I sample were all heavily ab¬ 
sorbed, this is not universally true of all 5 com¬ 
ponents; Hen 3-1591, EG And and 4 Dra are all 
lightly absorbed and detected below 1 keV. 

When compared with earlier data, all sources 
show changes in their intrinsic flux and N^f. In 
Paper I we proposed a scenario where the 5-type 
X-ray emission in symbiotics originates in the ac¬ 
cretion disk boundary layer, implying that long 
term changes in the X-ray flux are mostly due to 
changes in the accretion rate to the disk while 
changes in the soft X-rays can also be caused 
by variations in the amount of absorbing mate¬ 
rial. It is still unknown where the absorber is 
located. The days-to-week time scale changes 
i n Nff observed in t he 5A vpe prototype RT Cru 
( Luna et al. l2010l 12007 1 suggest that the ab¬ 
sorber is located close to the WD. However it is 
unclear if these changes are related and if so, how 
they are related, with the amount of material en¬ 
tering the accretion disk. In our study we witness 
that although high flux states encompass high Nh 
in the case of CD -28-3719 and 4 Dra, in the case 
of Hen 3-461, Njj is higher now while Ejc is lower 
than when observed with Swift (see Table 2 in Pa¬ 
per I). 

Our crude spectral fit indicates that the amount 
of absorption of the hard X-ray component in 
Hen 3-461 has increased considerably and possibly 
a separate soft spectral component has emerged. 
The intrinsic luminosity has decreased by about 
30% since the Swift observation in 2010. Al¬ 
though not contemporaneous UV data are avail¬ 
able, GALEX (NUV) observations taken one 
year before our Suzaku data indicate that the 
F( 7 V"= 1 - 86 x 10 “^^ erg cm“^ s“^, then Fuv/^x ^ 
0.7 (we quote the lower limit as the reddening 
for this source is unknown) and thus the accre¬ 
tion disk boundary layer seems to be still in the 
optically thin regime. In Paper I we proposed a 
scenario where the soft emission of /3/5-type ob¬ 
jects could be related with a colliding-wind region 
or jets. We could be witnessing the clear-up of a 
colliding-wind region. The equivalent width of the 
(unresolved) Fe Ka regio n of around 400 eVresem- 
bles the values found by Mukai et al. I ( 2007 1 on 
the well-known jet source with a two-component 
X-ray spectra, CH Cyg. 

Our best-fit model for Hen 3-1591 points out 
to the presence of a white dwarf as the accreting 


























Residuals normalized counts S"' keV' 


Table 3 

Timing analysis 


Object Bin size [s] Sexp s Ratio Sexp/average s/average 


4 Dra(4050) . . . 

160 

0.05 

0.17 

3.32 

15 

50 

4 Dra(4060) . . . 

160 

0.08 

0.39 

4.72 

8 

38 

Hen 3-1591 . . . 

160 

0.05 

0.05 

1.02 

52 

53 

4 Dra(4050) . . . 

1600 

0.02 

0.14 

7.27 

6 

43 

4 Dra(4060) . . . 

1600 

0.03 

0.33 

10.26 

3 

32 

CD-28 3719 . . . 

1600 

0.01 

0.01 

1.08 

47 

50 

Hen 3-1591 . . . 

1600 

0.02 

0.03 

1.30 

20 

26 

4 Dra(4050) . . . 

3680 

0.01 

0.14 

9.05 

4 

40 

4 Dra(4060) . . . 

3680 

0.02 

0.32 

13.71 

2 

30 

CD-28 3719 . . . 

3680 

0.01 

0.01 

1.24 

39 

48 

EG And. 

3680 

0.01 

0.01 

1.18 

52 

62 


®'We list only those bin sizes for which s > Sexp 



Fig. 2.— TS'C'T/GIS spectra of Hen 3-1591. The 
full line shows the best-fit model described in Sec¬ 
tion [31 


compact object instead of a neutron star, i.e. a 
white dwarf symbiotic instead of a symbiotic X- 
ray binary as previously though. Evidence sup¬ 
porting this scenario also comes from the opti¬ 
cal spectra which shows many emission lines of 
a variety of ions, similar to that of a planetary 
nebula, where the strong UV radiation field from 
the WD photoionize s the surrounding nebulae. 
Hedrick fc Sokoloski ( 2004l l detected the flicker¬ 


ing behavior in the B band data of this object, 
which suggests that the blue light is from an ac¬ 
cretion disk. Hen 3-1591 belongs to a rare subclass 
of d'-type yellow symbiotics, with a dusty IR con¬ 
tinuum and they belong to the Galactic disk, so 
low metallicity is not implicated. They are inter¬ 
preted as systems in which the hot component has 
recently evolved from the AGB to the white dwarf 
stage. In this interpretation, the dust is from the 
mass lost by the AGB star, and the nebulosity 
is in fact the planetary nebula (PNe), and nei¬ 
ther is due to the present-day giant (see Jorissen 
et al. 2005 and references therein). If we com¬ 
are with the well-known classical novae GK Per 


I 


Bode et al. Ill987rl and V458 Vul ( Wesson et al 


I 2 OO 8 I I. both claimed to be a classical nova inside a 
planetary nebula and if the PNe phase is as short 
as thought (10-20 kvr: lBadenes et aL~ ( 2015lH . we 
might question the white dwarf nature of the ac- 
cretor which should be still a hot subdwarf. 
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Hen 3-1591 also shows the barium syndrome, 
i.e. overabundances of s-process elements and the 
presence of single-ionized barium which cannot be 
explaine d if the red giant is not p art of a binary 


syste m ( Jorissen fc Mayor 1992t Jorissen et al. 


200511 . The now-observed red giant had its pho¬ 


tosphere polluted by s-processed elements by the 
previously-AGB companion, that now should be 
a WD. Thus, our high signal-to-noise, broadband 
Suzaku spectrum adds more support to the pres¬ 
ence of a WD or hot subdwarf. If we use the 
plasma temperature from the cooling flow model 
to derive the WD mass, we obtain MvyD= 0.45 
Mq, on the low-side of the WD mass distribution. 


The high temperature of the plasma strongly 
suggest that in Hen 3-1591 the X-ray emission 
arise in an accretion disk boundary layer instead 
of a colliding wind region. The strong-shock con¬ 
dition implies wind speed’s of around 3,000 km 
s“^ for the observed temperatures, and such high 
speeds outflows or winds had not been detected 
so far in Hen 3-1591 or any other symbiotic. The 
lack of UV data does not allow to use the ratio 
of UV and X-ray fluxes as a proxy for the opti¬ 
cal depth of the boundary layer. The decrease of 
temperature and luminosity (i.e. accretion rate) 
between ASCA and Suzaku observations, how¬ 
ever, suggest that the optical depth of the X-ray 
emitting plasma has changed, being higher dur¬ 
ing Suzaku observation and thus we might have 
observed an smaller optically thin portion of the 
boundary layer. 


Based on the high temperature found for the 
X-ray emitting plasma, we can also rule out col¬ 
liding winds as the source of X-rays in EG And. 
This s ource was identified as a /3-type symbi¬ 
otic bv iMiierset et al. (1997). Those authors de¬ 
rived a temperature for the X-ray emitting plasma 
1.3±0.5 keV using ROSAT PSPG data. This value 
is much lower than the temperature found from 
our modeling of the Suzaku spectrum (5-10 keV) 
and is much likely due to the lack of sensitivity 
above 2 keV of the ROSAT spectrum. The ve¬ 
locity implied by the best fit to the Suzaku spec¬ 


trum is around a few thousand km s 


This is 


much faster than the highest velocity line features 
(^700 km s“^, G IV 1548, 1550Aabsorption fea- 
tu res) found in UV spect ra from FUSE and STIS 


bv ICrowlev et al. ' ( 2008ll . and difficult to explain 


in a low mass white dwarf like the one in EG And 


{Mwd -0.4 M^: lKolb et al. 11200411 . 


Instead, we propose that the X-rays originate 
in the accretion disk boundary layer around a 
low mass WD. We estimate the WD mass tak¬ 
ing into account the relationship between the 
kT from the spectral modeling {kT=7t.\ keV), 
kT—(l/10)/rmp(GMwD/Rivn), where is the 
mean molecular weight and rrip is the proton 
mass and the WD m ass-ra dius relationship from 
Hansen fc KawaleTI ( 1994ll . we obtain a mass 


Mwd —0.5 Mq. 

If the hot-component lu minosities rep o rted i n 


the 


Kolb et al. 


literatur e (16 -400 Lq I Vo gel et al. (|l992ll ; 
(l2004lU are due to an accretion- 


powered white dwarf, then the implied accretion 
rates are in the range 10 “® to 10 “^ Mq yr“^ for a 
0.4 M 0 WD. This is squarely in the regime where 
the boundary layer is exp ected to be optically 


thick according to models bv lPooham fc Naravan 
I I 995 II . The X-ray luminosity is lower than L/iot, 
suggesting that the X-rays are produced in a re¬ 
gion that remains optically thin at the outer sur¬ 
face of the mostly optically thick boundary layer. 
Although we expect X-ray emission due to accre¬ 
tion to be highly variable, EG And is very faint 
and we do not detect enough photons to be sensi¬ 
tive to low-level variability {sexp/average — 50%). 

Our Suzaku observations of 4 Dra, when consid¬ 
ered along with previous UV observations ( Skopall 
l2005ll . indicate that 4 Dra is an accretion-powered 
symbiotic (i.e., no quasi-steady shell burning on 
the surface of the WD) in which mass transfer 
at —10“® Mq yr“^ produces an optically thick 
boundary layer around a —0.6 Mq WD. We ob¬ 
tain this rough estimate of Mwd from our infer- 
ence that the boundary la yer is optically thick. 
Popham fc Naravan ( 1995ll suggest that the ac¬ 


cretion rates needed to produce a predominantly 
optically thick boundary layer around WDs of 
masses 1.0, 0.8, and 0.6 Mq are greater than 10“^, 
a few times 10“®, and — 10“® M 0 yr“^, respec¬ 
tively. Taking the UV flux at the time of the 
Suzaku observation to be comparable to the —10 
L 0 detected by lUE, the BL being optically thick 
suggests that the WD mass is less than about 0.6 
M 0 (for which the lUE luminosity corresponds to 
accretion rates of —10“® M 0 yr“^). In contrast 
to higher-mass WDs with optically thick bound¬ 
ary layers (such as RS Oph and SS Gyg), the small 
amount of X-ray emitting plasma associated with 


10 























































the optically thick boundary layer around 4 Dra is 
just as hot as would be expected for an optically 
thin boundary layer. 


We emphasize that our conclusions about the 
optical depth of the accretion disk boundary 
layer and its implications on the WD mass de¬ 
pend on the assumed viscosity parameter, a and 
other assumptions on the models described by 


Popham fc Naravan (1199511 . As noted in Paper 


I, a 30% change in a leads to a factor-of-a-few 
change on the threshold accretion rate. Allowing 
for this, the white dwarf in 4 Dra is smaller than 
about 0.7 Mq, still a low mass. 


Both fitting approaches discussed in Section 
13.1.11 yield equally possible scenarios. In the case 
that both spectra are modeled independently, we 
can understand that the measured temperatures 
are different because the optical depth of the emit¬ 
ting region has changed, i.e. during the earliest 
observation the boundary layer was thicker and 
the measured temperature and fluxes of the X- 
ray emitting plasma were lower. On the other 
hand, if we assume that during both observations 
the plasma temperature was the same (the optical 
depth of the boundary layer and the mass of the 
WD have not changed between observations) we 
still obtain acceptable fits with an increase in the 
amount of absorbing material and the intrinsic lu¬ 
minosity during the last observation, perhaps due 
to an increase in the accretion rate. 


The orbital soluti o n for 4 Dra was stud¬ 
ied by IPamaev et al.l ( 20091) . conducting a ra¬ 
dial velocity study of the M giant and found a 
1,703±3 days period, an eccentric orbit (e=0.3) 
and To[MJD]= 53204±19. The ROSAT observa¬ 
tions thus occurred at orbital phases 0=0.33, 0.15 
and 0.61, while S'wzafcws ’ occurred at 0^°®°=O.23 
and 0"‘°®°=O.57. If the X-ray absorbing material 
is tied to the orbital motion of the WD, thus the 
absorption columns obtained from the fits of the 
Suzaku observations should not be very different 
given than during ObsID 4050 the system was near 
quadrature while during ObsID 4060 the system 
was in inferior conjunction, thus from the values 
in Tableland since an undisturbed M giant wind 
would account for only ^^ a few times 10^^ cm' 


( van den Berg et al. 20061) . we can conclude tha'; 


the absorber material is located relatively near the 
white dwarf and probably intimately connected 
with the physics of accretion. 
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